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Talbot effect in the space-time evolution of matter waves is analyzed 
and shown that the matter waves at relativistic and non-relativistic veloc- 



ities exhibit coherence beyond the grating and display Talbot self-imaging. 
The grating is realized by considering equidistant narrowly peaked coherent 
^ ■ probability distributions at the onset. Theoretical frame work is based on a 

quantum particle at relativistic and non-relativistic velocities in a one dimen- 
sional box. The matter waves, in their evolution, display Talbot self-imaging 
and fractional self-imaging as restructuring, respectively, at Talbot length 
and fractional Talbot lengths. 
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I. Introduction 

The Talbot Effect in classical optics, observed in 1836 by Henry Fox Tal- 
bot, is a near-field diffraction effect [1]. As a laterally periodic coherent 
wave is incident on a diffraction grating, the grating image repeats at regular 
distances away from the grating plane. The regular distance is called the 
Talbot Length. At fractions of the Talbot length we find fractions of the 
grating images [2]. The phenomenon has applications in optical metrology, 
photolithography, optical computers [3], as well as, in data compression [4] 
electron optics and microscopy [5]. Sir Michael Berry and Susanne Klein 
have shown a relation of Talbot effect with quantum revival phenomenon of 
a particle [6] . Later, in another seminal paper a close relavance between Tal- 
bot effect with light and quantum revivals of atoms in a one dimensional box 
was developed [7]. Integer and fractional nonlinear Talbot effect is recently 
reported in nonlinear photonic crystals [8]. In atom optics the effect is seen 
using atoms with non-relativistic velocities, displaying self imaging at Talbot 
lengths [9]. Particle in a box has been treated as playground to study quan- 
tum revivals [10, 11] and space-time dynamics which weaves quantum carpets 
as a result of delicate interference between eigen modes of a system [12, 13]. 
For a slightly relativistic particle in a one dimensional box, quantum carpets 
are explained earlier by using Green's function where slight relativistic cor- 
rections appear in times of revivals [14]. One dimensional Dirac equation is 
used to treat a relativistic particle along a circle [15]. In this contribution, 





Figure 1: Left panel shows that matter waves passing a few-slit grating with 
non-relativistic velocity, v, display interference pattern but no Talbot self- 
imaging. Right panel describes that matter wave passing a diffraction grating, 
defined mathematically by equi-distant narrowly-peaked probability distribu- 
tions, display Talbot fringes even in early evolution. In the early evolution 
individual distributions after respective initial quantum dispersion in free 
space interfere with each other. Here, grating spacing s = Lj 10, grating 
elements widths A = 0.001L, z = vt, where L is the length of the grating 
and t is evolution time. 



I explain Talbot effect for matter waves at non-relativistic and relativistic 
velocities. Theoretical treatment is conceptually based on the evolution of 
a quantum particle in a one dimensional box, which, serves as a finite uni- 
verse model to explain Talbot effect. I show that the matter waves manifest 
interference leading to restructuring at times which determine self-imaging 
of Talbot. Furthermore, I illustrate that matter waves with relativistic ve- 
locities display unique Talbot carpets in space-time evolution, different from 
classical optics. Hence, I provide a rigorous treatment to the problem based 
on relativistic quantum mechanics, first, to understand Talbot effect with 
matter waves and, second, to explain the evolution of a relativistic particle 
in a one dimensional box. 

In order to explain Talbot self imaging for matter waves, let us consider an 
initial array of equidistant narrowly peaked probability densities in a ID box, 



which serves as equidistant wave-fronts emerging from grating. In order to 
observe Talbot effect, the distance between two wave-fronts, which is suitably 
regarded as grating-distance, is taken much smaller than the length of the 
box. During its propagation the array displays quantum interference, for a 
very early evolution, see Fig. 1 (right panel). The probability density for the 
matter wave at position, x, and time, t, is calculated as 

W(x,t)= i dx' f dx"K(x,t\x',x")uUx')uj(x") (1) 
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and, therefore, follows from the integration of the kernel 

K(x,t\x',x") = G\x,t\x')G(x,t\x"). (2) 

Here, u = u>(x,t = 0) denotes the initial wave function and, G expresses 
Green's function, defined as, 

oo 

G(x, t\x') = J2 Vm{x) cp m (x')e- iE ^ h , (3) 

m=l 

written in terms of the energy eigen functions (p m and energy eigen val- 
ues, E m . The evolution of a particle of mass M in one dimensional box of 
length L, along x— axis, may be expressed by the one-dimensional relativis- 



tic Hamiltonian, H = J(Mc 2 ) 2 + (p x c) 2 — Mc 2 . In the limit c — > oo the 

Hamiltonian H reduces to non-relativistic form, whereas, retaining relativis- 

tic corrections to order 1/c 2 provides slightly relativistic Hamiltonian. The 

solution, cp(p, x) = yVsin( ,§/h), satisfies the equation < H — E p \ cp(p, x) = 0, 

v 
in non-relativistic and slightly relativistic limits [14], where p = hk = rrnih/L, 

and M is normalization constant. 



II. Talbot effect with relativistic matter waves 

In relativistic case we solve the Dirac equation for one dimensional box, that 

is, 

H<p = (a z p x c + /3Mc 2 ))y? = E<p, (4) 

where tp denotes two-component spinor depending upon x, belongs to [0, L] [16], 
and a z and (3 are the matrices 

az= {l a o) ' ^ = (o 7 -/)- 

Here, a z is Pauli matrix and / is unit matrix. The solution of the Dirac 



equation within the box is 



2cos ™p_ I ) x 



where, the constant M is obtained from the normalization condition. In 
addition, we write, 

6m = arctan (i^%) Pm = ETW^- (6) 

Moreover, x is an arbitrary two-component normalized spinor, i.e. x^X = 1> 
and 



E m = Mc^l + (2mg) 2 . (7) 

The controlling parameter q is defined as, q = (fff ) /4L = A C /4L = t>/c, 
a ratio between Compton wavelength A c and four times the length of the 



box, or a ratio between effective velocity v and speed of light, c. In order to 
understand the space-time dynamics, we write the Green's function for the 
component which follows sinusoidal dependence on position, x, that leads to 
the calculation of the kernel, viz., 



K = D r 



D r 



D r 



D 



r\ i i Ji 



(8) 



Here, we use scaled coordinates £ = x/L and r = t/T, in addition £„ 
TE m /nH and T = 4ML 2 /vrft, and 

oo 

D r (+,+) = D r (+r ] ,+C)= J2 e^ m "-<) /+,+(£) 

m,n=— oo 
iirl (m-n)^-(TMc 2 /-Kh){yJl+(2mq) 2 -^l+{2nq) 2 }T\ 



(9) 



Here, /+,+ = 1, however, /__ = e i(.^(m)S m - S gn(n)5 n ) ^ j- + _ = e -is g n(n)«5 n? and 
/_ j+ = e lsgn(m ) <5m . In exponent sgn defines signum function. 
The double series transform 
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m,n=— co 



j,l=-oo 



2 0' + P) ' 2 ( J ~ p ) 
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provides the following expression for _D r (rj, (), i.e., 



Dr(+,+)= E El- 1 )'' 

Z=— oo j=l 



3 J7TJ 



1-C 



e -"j 2 



■■ T|-C 



X+ + e-^'^Z- + E x < 



-o/, 



J=— oo 



where 



(10) 
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I± — dpcos (nipp ± a{s + (p) — S-(p)}) 
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1< 



(i/ 



dpcos(nipp). 



In order to simplify the expressions here we use the following abbreviations: 



i> 



v+( 



+ £ + I; a = Mc 2 T/irh; and s±(p) = Jl + q 2 (j ± p) 2 . 



I solve the integrals, Z±, and substitute the expressions in Eq. (1). Hence, 
for a matter wave with finite initial distribution, A, much smaller than the 
grating spacing, s, the probability density becomes, 



oc 



oo oo 



w(t,r)=Ar*^= E E 



en 



-oo 



~| 6 7r e A« 



t=-°° |=oo 



n,m=~oo ki,tt2=—k 



,in +r .(+) 



+ (-i: 



- t 7+(fc 1 ,A; 2 )(e tfl -ri- ) + (-l) 



i+m — iQ_i_T-i(+) 



n+m — i 



r_ 



fi-p(-) 



(12) 



Here, I introduced J-{k\, k 2 ) = J m (7_)e OTr - 7— a - and J+{h\, k 2 ) = </ m (7+) e ~" r: '^~ 
Furthermore, /ci and /c 2 are integers and /c corresponds to total number of 
initial narrowly peaked probability densities at grating spacing, s [17]. Fur- 
thermore, Jl± = J7r(£ ± 2r/g + £ ± (k± + fc 2 )/2), and r± is defined as 



7T J-oo 



(13) 



(±) 



The parameter a^ = n(ki + k 2 + 2ijnA 2 )/2q and Xo = ^(C =t r /<7 + ^ + 
kl 1 2 k2 )/q. The four terms in the expression for K in Eq. (8) pick four different 
corresponding combination of signs for ki and fc 2 . Hence, the full expression 
becomes four-fold to what we have in Eq. (12). 

The final expression though difficult in its structure is easy to analyze. 
As discussed above n and m satisfy the equation n + m + 2 > 0, keeping 
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in view minimum criteria we choose n = —2 and m — 0. Furthermore, we 
neglect imaginary part in a as A 2 << 1, that simplifies the expression for 
W(£, r), as it becomes 



oo °° °° 



|=— oo |=oo 



k 
X 

k\=—k 



]T [(J (7 + )cosQ_r(-)H-Jo( 7 -)cosQ + r( + )]. (14) 



For simplicity I consider k\ = k^. Here Jo(7±) is a function of time, and 
cosf2± defines the space-time structures. Furthermore, solving Eq. (13) in 
the vicinity of Xo , we get, 

r(±) ( ± ) 

Ux^) - ^{in-i{xt } ) + 4 + i(xi ±} )), (15) 



x 2ttA 
which is independent of j . 

With the help of asymptotic expression Jq(z) ~ J2/nzcosz, we get 

Jo( 7 ±) cosfi T = cos(j7r(e =F (2/g =F 1/0V))t + 0)- (16) 

Note that for all the values of j we have periodic behavior both in space 
and in time. In (£, r) plane there exists maximum probability lines with 
slope =F(2/g =p l/(j7rg 2 )). It is interesting to note that the j dependent term 
in the expression for slope becomes negligibly small for q > 1, however, for 
smaller q values it contributes, removes degeneracy, and support many canals 
and ridges behavior corresponding to different j values. Hence, for q — 1 
and higher, we effectively find higher probabilities following straight lines 
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with slope ±g/2 for all values of j, which manifests quantum coherence at 
relativistic scale. The analytical results show a good comparison with exact 
numerical results. For the component of Dirac solution, given in Eq. (5), 
which does not satisfy the boundary conditions we get similar results for 
probability density with a difference of phases between constituent terms 
corresponding to D r 's, given in Eq. (8). 

III. Talbot effect with non-relativistic matter 
waves 

As the particle evolves with non-relativistic velocity, the corresponding ex- 
pression for D r (+,+) becomes 

£ (-iy l e^^5(l + Z-2 3 T+^±^\. (17) 

j,l=—oo V / 

The probability density, as calculated following Eq. (1), reveals 

OO ±00 g £ g 

$=— oo |=±oo 
ki,k2=—k 

+ ^ E e-^(^ + ^). (18) 

l=— oo ki,k2=—k 

As in relativistic case, the four terms in Eq. (8), respectively, take one corre- 
sponding sign combination for hi and k 2 . Hence, the complete expression for 
W is composed of four similar terms, given in Eq. (18), with four different 
combination of signs for k\ and k 2 . Beyond r = 0, a jet of sub-distributions 



emerges following straight lines of slope 2j in (£, r) plane and shifted in time 
by I ± fcl ^ fc2 . We note that the dominant values of j which contribute in 
defining the probability distribution as a function of space and time have 
width, l/y^vrA. For the reason smaller value of A corresponds to more 
straight lines which define richer space time structures. Hence we find, for 
each straight line with slope 2j there exists a counterpart with slope — 2j and 
their interference carves carpet structures, with a periodicity in time. The 
last summation term in Eq. (18), explicitly written for j = 0, corresponds to 
straight lines with zero slope, along the r axis, for / and ± kl ± k2 , 

Matter waves with non-relativistic velocity display Talbot self-imaging, 
similar to light waves in classical optics, shown in Fig. 2 (upper panel). The 
space time dynamics for non-relativistic matter waves is defined by Eq. (18), 
which depicts that the quantum Talbot carpets are more dense in structure 
as compared with quantum carpets for a particle in a box, as there exist 
4k(k + 1) more low and high probability density lines in quantum Talbot 
carpets, respectively, named as canals and ridges. Here, k is the total number 
of initial probability densities or wave-fronts which define grating. The self- 
imaging takes place at Talbot length, zt = 2A^s 2 [2], and fractional imaging 
occurs at fraction of the distance. Here, s is the grating's spacing and X^b is 
the de-Broglie wavelength of the matter waves, defined as XdB = 2irh/Mv = 
(Tg^a) -1 , as the particle moves with effective velocity v, and T is the time 
of quantum revival. 





Figure 2: Talbot self-imaging for matter waves with non-relativistic velocity 
(upper panel). At fraction of Talbot length we note fraction of grating im- 
ages. For matter waves with relativistic velocity (lower panel) we find simpler 
Talbot carpet. Since for each initial wave-let there appears two symmetric 
non-zero probability density lines, explained in the text, the final Talbot car- 
pet for matter wave is different from non-relativistic one and, that for e. m. 
waves in classical optics. The Talbot length, however, is larger now as com- 
pared with non-relativistic case due to time dilation effect. The parametric 
values are the same as in Fig. 1. 
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Matter waves with relativistic velocities display self imaging in complete 
contrast to non-relativistic case, and to e.m. waves in classical optics. As dis- 
cussed above, Eq. (14), in relativistic case there are only two possible non-zero 
probability density lines in (£, r) plane, and wave packet intersect each other 
periodically, thus, defining full restructuring of original wavelets. This leads 
to self-imaging of the grating at Talbot lengths. Whereas at any intermediate 
time there occurs self-imaging twice of the grating structure, see Fig. 2(lower 
panel). We defines the Talbot length, z T = 2A^s 2 = 2(T re/ ^)- 1 . Hence, in 
relativistic case Talbot carpet is simpler in structure, as compared with clas- 
sical optics, however Talbot length is longer due to time dilation relativistic 
correction [18]. 

IV. Conclusions 

I conclude that in the present contribution, I develop the theory of Talbot 
effect for matter waves with relativistic and non-relativistic velocities. I show 
that quantum coherence prevails and drastically modifies the interference be- 
havior in relativistic case which results in an exact Talbot imaging, available 
in non-relativistic quantum mechanics as well, however, following different 
Talbot carpets. The mathematical analysis provides theoretical understand- 
ing of a relativistic particle in one dimensional box as well. Keeping in view 
latest technological development, the analytical results can be realized using 
cold neutral atoms [9], and in recent experiments with electrons [19]. 
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